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Abstract 

We report on a nontrivial bosonization scheme for spin operators. It is shown that in the 
large N limit, at infinite temperature, the operators Y2k=i ^fci/V^V behave like the creation and 
annihilation operators, a' and a, corresponding to a harmonic oscillator in thermal equilibrium, 
whose temperature and frequency are related by huj/ksT = In 3. The z component is found to be 
equivalent to the position variable of another harmonic oscillator occupying its ground Gaussian 
state at zero temperature. The obtained results are applied to the Heisenberg XY Hamiltonian at 
finite temperature. 
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I. INTRODUCTION 



The interest in bosonic representations of spin operators arises from the necessity of 
finding alternative practical tools to the diagrammatic technique for spins which turns out 
to be quite complicated [Li]. By using bosonic operators, one can, for instance, evaluate the 
functional integrals and afterwards calculate the different contributions diagrammatically (2)- 
8. Moreover, spin-wave theory is based on the possibility of transforming spin operators into 
new bosonic operators, which are much more easier to deal with [5]. This theory has known 
great success in the study of magnetic properties of many materials at low temperatures; 
it led, in particular, to the notion of magnon, the elementary collective excitation of spins 
in the crystal lattice. Many important spin Hamiltonians, such as the Heisenberg model 
Hamiltonian, can be exactly diagonalized at low excitations with respect to the bosonic 
degrees of freedom using certain transformations. 

Among the most used transformations in spin-wave theory we cite the Dyson Ql and the 
Holstein-Primakoff transformations. The latter maps spin-| operators as follows [7|: 

s + = a 1 \/ \ — a^a, (1) 
s_ = vT— a^aa, (2) 
s z = -- + at a. (3) 

Here and throughout the text, unless otherwise stated, we take h=l. For a set of N spins, 
the scaled total raising operator is thus given by 

N 



y/N VN 



LZXfcV 1 - a W- ( 4 ) 



fc=i 

.t. 



At low temperatures and low excitations, i.e. when (a k ak) << 1> we ma y ; as a first approx- 
imation, write 

S+ * N 



^£4- (5) 



For large N, the Fourier transforms of the bosonic operators and a\ are given by 
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These are also bosonic operators, since they verify the commutation relation [ai, a\,] = Set'- 
We can thus conclude that at low excitations (temperatures), and large N, 



which is indeed a beautiful result. However, at higher excitations, the square root in equa- 
tions ([I])-© should be expanded in powers of ala, depending on the order of the adopted 
approximation. This gives, for instance, 

7^ = 7^h 2 s r * 

In this paper we address the problem of finding a suitable bosonic representation similar 
to ([7j) that is valid at high temperatures. As we shall see, the cost we have to pay for 
preserving the above form of the transformation resides in the introduction of effective 
temperatures that are generally different from the actual ones, together with a deformation 
of the functions of interest. 

The paper is organized as follows. In Section [Til we prove some fundamental results in 
connection with the addition of spin operators. The focus is on the trace properties of the 
suitably scaled operators to ensure good statistical behaviour in the large N limit. Our main 
result is presented in Section [TTT1 where we explicitly establish the exact form of the mapping 
between the spin and bosonic operators. The paper is ended with a brief discussion. 

II. PRELIMINARY RESULTS 

Let us begin by proving some results about the behaviour of the spin operators S a / VN = 
Ylf=i & ia/ "^A (a = x,y,z) as A — > oo js _ 12]. By the multinomial theorem we have that 



tr(s a y= - , r j~ - , tr ^ ® ^ ® • ■ ■ ® § N N j(t - E r • ( 9 ) 

rm-TN 1 ' 2 ' N ' k 
Note that we are using the Kronecker delta symbol. We have to distinguish between two 

possible cases, depending on whether i is odd or even. In the former case, we see that there 

should always be an operator, in each term of the right-hand side of the above equation, 

appearing with an odd power, say, 2n + 1 where n G Z+. But since trs^ re+1 = for all n, we 

conclude that 

/ £ \ 2£+l 

tr(4M =0, WGZ+. (10) 



When the power of the operator is even, 

ri^—rjv A;=l ^ k ' 

One can easily see that the sum in the right-hand side of equation ffTTl) yields a polynomial in 
iV of degree £, the main contribution of which comes from the term in which all the operators 
Sk a contribute equally to the product. This happens when all the operators appear with the 
same power; all the other possibilities yield polynomials of degree less than t. But due to 
the restriction imposed on the sum of the powers of the operators, we conclude that: 

2/1 / e N \ 

tr ^ = 2 ! X2! ... 2! Eii^ n tri +^-i(ao, (12) 

J > n e \k=i k=e+i / u t [i,2,-N] 

l terms 

where the products are evaluated for all the possible partitions ILj[l, 2, • • • , N] of N elements 
into subsets of I elements, and Qi-i(N) is a polynomial in N of degree at most equal to I — 1. 
Next we remark that for spin- 1/2 operators, s 2 a = meaning that trs^ = ~, independent 
of the partition He. Hence we can write 



tr( . j2< = ^ N(N . mN ^... { N- t+l) {tisl)t{tih ^ + Qt _ i{N) 



2 N N e 




fe=l 



(13) 



lirn 2- N tr(S a /VN) 2e = , a ee x, y, z. (14) 



It follows that 



Now we shall prove that as iV — > oo the operators S a /y N and Sp/y N become uncorrelated 
(i.e. independent) with respect to the tracial state when a ^ (3. We have 

m\l\ 



, . , r 1 \r 2 \---r N \k 1 \k 2 \---k N \ 

k 1 ,k 2 ,---k N ,r 1 ,r 2 ---r N 



, ri\r 2 \- ■ -r N \ki\k 2 \- ■ -k N \ 

fcl,fc2,-"fc]V,»*l,»*2"-rjv 

x tr( s? a . s%) ® (s* 2 l.s%) ® . . . ® 
x 5 ( m 



(m-^r^5^-^^. (15) 



Once again, if one of the powers is odd, then in each term in the right-hand side of the latter 
equation there should be an operator §j appearing with an odd power; this means that its 
trace is zero, and hence the overall trace vanishes as well. This gives 

ti(S 2 a e+1 Sf) = ti(S a ) 2e+1 ti(S(s) m = 0, Vme Z+. (16) 

Consider now the case in which the powers of the total spin operators are even and make 
the substitution m — > 2m, £ — > 2£ in equation f|T5|) . It follows that the right-hand side of 
the latter equation yields a polynomial in iV of degree £ + m; we shall assume that £ > m. 
In this instance, the leading contribution to the sum comes from the terms in which all the 
spins equally contribute to the trace, which occurs only if all their powers are equal to 2, 
taking, obviously, into account all the possible partitions of N elements into disjoint subsets 
of £ and m elements. Since we have assumed that £ > m, we should have a contribution to 
the trace of m operators in the form s 2 a s 2 q, and £ — m in the form s%; the remaining N — £ 
traces come from the identity matrix. More precisely: 

(2m)\(2£)\ /m 



£ terms m terms 

+ P e+m -x(N), (17) 



where Pe +m _i(N) is a polynomial in iV of degree at most equal to £ + m — 1. Taking into 
account the fundamental properties of spin-^ operators, we find that 



tr 



//dy^/dyKl (2m)!(2l)! N(N - 1)(N - 2) ■ ■ ■ (N - m + 1) 
N(N-l)(N-2)---N(N-e+l) 2 N ~" 

x j\ x 8 m 2 e-m + n+m-i{Jy) 



N 



l+m 



«> ! n (l -rri (l -!) +0 (i 

k=l k'=l 



2 3e 2 3m £\m\ A1 V A^yilV NJ \N 

k=l ' ■ ■ 

As a direct result of equation ffTSl) we deduce that 



;is) 



Ita 2-tr {(SjVN?^,VNr} = MB, (19) 



which means that 



lim 2-*tr \(S Q /VN) 2m (Sp/VN) 2e \ = lim 2 _Jv tr (^/v/iV) 2 ™ lim 2~ N tr (S^/Vn) 2£ . 

(20) 



Equation ()14p determines the explicit form of the moments of the random variable, i], 
associated with the spectral decomposition of the operators S a / \/N as N tends to infinity. 
To find the expression of the corresponding probability density function p(rj), we first note 
that the characteristic function is given in terms of the moments (r] n ) by 

00 / n\ 

so that 



n=0 



n=0 v ' n=0 

This yields 

= e"' 2/8 . (23) 

The probability density function is nothing but the Fourier transform of the characteristic 
function: 

1 r°° 

p(rj) = — / <&(t)e- ltri dt. (24) 



P(V) = \l -e~ 2r , (25) 



By direct substitution of the expression of $(£) into the latter equation, we find that 

1 

7T 

that is, a Gaussian probability distribution with zero mean and variance a = ~. 
A corollary of the above result is that for any well-behaved functions /, g: 



lim 2T N trf(-%) = Jl ^ f{ V )e-*dn, (26) 

* * +00 

lim 2-"tr0(4f=, -%) = -[[ g(v,n)e- 2 ^ + ^dr)d», a, = x,y,z. (27) 

—00 

If, instead of the operators 5^ and S y , we use the lowering and raising operators S+ and S-, 
which are related to each other by 

& = i(£ + + £_), = (28) 

the integral in equation (|27|) should be expressed in terms of the random variables z* and £ 
corresponding to S + and SL, respectively One can easily verify that 

lim 2~ N trg( 4=) = - / I g(z* , z)e~ 2 ^ dz*dz. (29) 

This being explicitly proved, we are ready now to formulate our main result in the next 
section. 
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III. BOSONIZATION 



The crucial observation is that equation (|29|) is nothing but the coherent state repre- 
sentation of the mean value of the function operator g written in normal order (which we 
denote by Afg): 

-J J^g{z\z)e- 2 ^ 2 dz*dz = J P{z* , z)Mg{z* , z)dz*dz. (30) 

Here P(z*, z) is the P-representation 13J of the density matrix p of an equivalent bosonic sys- 
tem (harmonic oscillator, as we shall see bellow) whose creation and annihilation operators 
are denoted by a) and a, respectively. Explicitly, we have 

P( z * ? z ) = J2(n\pS{z* - a?)5(z - a)\n) (31) 

n 

which, obviously, is normalized to unity: 

J J P(z*,z)dz*dz = 1. (32) 

J z* J z 

It follows that the density matrix p is given by 

p= I P(z*,z)\z)(z\dz*dz, (33) 



with the coherent states 



oo r 

z r ' 



Whence 



ml 

m=0 



or r 7 n 7 *m 

(n\p\m) = - dz* dze- 3 ^^== (35) 
7T J z * J z vnlml 

r) i>2ir poo n+m 

= - d</> e~ 3r2 - e^-^rdr (36) 

7T Jo Jo Vnlml 
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(37) 



"""3n+l ' 

where use has been made of the polar coordinates (r,<f>). In fact we can rewrite the above 
result in the from 

Wl"> = (1 +i)n + r ( 38 ) 
This reminds us the form of the occupation number representation of the density matrix of 
a harmonic oscillator in thermal equilibrium at temperature T, namely, 

P = n nK\ -li \ n ) ( n \ ( 39 ) 

r L^i (l _|_ ( n ))n+i 1 /x 1 v ; 



where the mean value is given in terms of the natural frequency of the oscillator by 

1 



(n) = tr(pa^a) 



(DK B J 



Comparing equations fl38|) and (1391 . and using equation fHUj) . we obtain 

1 fyjj 



(n) 



2' k R T 



In 3, 



meaning that 



p = e 



-ln3(at a +i) 



IZ, Z = 4- 



Our main result may thus be formulated as follows (h = 1): 
Theorem 1 Given any well-behaved function f, we have: 



lim tis\ ^ 
.v— x I ^2*^ 

p s (T=oo) 



=tr< e 



la3 ^ a+ ^Mf{a^a)\ 



(40) 



(41) 



(42) 



(43) 



We conclude that at infinite temperature and large N, the raising and lowering scaled spin 
operators S + /y/N and S_/y/N behave like the creation and annihilation operators a) and a 
of a harmonic oscillator in thermal equilibrium whose frequency and temperature are related 
by fiw/kBT = In 3. 



Worked example 

Let us illustrate the above result by considering the operator 



/ 



For N = 2000 we find that 



N VN 



\ = (S+S- + S-S + ) 5 /N 5/2 . 



I/V „ / S-L. S— 



119.670. 



On the other hand, 



(44) 



(45) 



Aff(a\ a) = 2V 5 a 5 = 32[(a f a) 5 - 10(a f a) 4 + 35(a f a) 3 - 50(a f a) 2 + 24a f a]. (46) 



Taking into account the fact that 

2 

V 71=0 



^ e -ln3(„ + i) n fe 



2 /l 
-Li_t. — 

3 V3 



(47) 
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where Li m (x) denotes the polylogarithmic function, we find that 

—= 32e- ln3(n+ ^(n 5 - 10n 4 + 35n 3 - 50n 2 + 2An) = 120 (48) 

V 3 n=0 

as should be. 

What about the z component? To determine the equivalent degree of freedom, let us 
have a look at the probability density 

x i-> p(x) = \ j^e~ 2x2 (49) 

and compare it with the Gaussian probability distribution of a harmonic oscillator [see 
equation ( l26l) ] 

,2 1 __r (x-(x)y 



x\-+\ib(x)\ = . exp<^ - , A / >, (50) 

1 /2tt(Ax) 2 I 2(Ax) 2 ( 



We find that 



"+00 

(/(*))= / f(x)mx)\ 2 dx. (51) 



(x) = 0, Ax = ^. (52) 



For the ground state (that is a minimum uncertainty state), the variance is given in terms 



of the mass m and the frequency to of the oscillator by Ax = y 2mcj (remember that we set 
h = 1). Hence 

muj = 2. (53) 

This means that, as iV — y 00, the operator S z /y/~N behaves like the position variable of a 
harmonic oscillator (different from the above one) in its Gaussian ground state (T = 0) whose 
variance is equal to 1/2. The full spin system is thus mapped onto a bipartite bosonic system, 
consisting of two independent harmonic oscillators, one of which is in thermal equilibrium, 
while the other one is at zero temperature. 

Let us apply the above results to the Heisenberg XY Hamiltonian (7 is the coupling 
constant): 

H = ^{S+S- + S-S+). (54) 



The thermal expectation value of the function operator f(S+/y/N, S-/y/N) is given by 



As iV — > oo, equation (143]) implies that 

= tr e- ln3 ( atQ +^Are"^ (ata+aat V(gt,a) 

where 7 is now given in units of h. Taking into account the fact that 

n 

M{a ] a + aa ] ) n = 2 n (a t ) n a n = 2 n ^ £"(a f a)* (57) 

£=0 

where B™ denotes Stirling's numbers of the first kind [14], we find that 

00 „ 1 n 

Ve -*(.w> = E (z|)» 1 EB?(aV . (58) 

n=0 B ' 1=0 

By using the following formula for Stirling's numbers of the first kind: 

OO 71 n 

( 1+i )" = EEV' (59) 

n=0 k=0 

we deduce that 

M e -T^T^*+^) = {l _ 2 1 /(k B T)f a = e M1 ~^ Wa . (60) 



From equation ( 1561) . we find that the partition function of the corresponding bosonic system 
is given by 

(61) 



-in 



( i- 27 A B r )( ata +^) 



which is the partition function for a single mode harmonic oscillator in thermal equilibrium. 
The following conditions should however be satisfied: 

1 > T^L 1 > — r^r- ( 62 ) 
k B T' k B T y 1 

In the case of ferromagnetic interactions, 7 < 0, we find that 

k B T > M- (63) 

This gives the lower bound on the temperature above which the bosonization is valid. To 
proceed further, let us apply equation ([7]) to the low excitation sector of the Hamiltonian; 
one finds 

H — T( a o Q; o + h.c) — 2j(ala + -). (64) 
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Hence, 



k B T> 



'0 



Hujq = 2\j\. 



(65) 



2 



We may thus define an effective temperature for the harmonic oscillator as follows: 

2\j\ 

^bTcS = - — j 3 r-- 

Finally, regarding the function /, it can be mapped according to 



IV. DISCUSSION AND CONCLUDING REMARKS 

In conclusion we have established a bosonic representation for scaled total spin operators 
at high temperatures. It should be pointed out that this bosonization scheme does not 
preserve the temperature, since the spin system, originally taken at infinite temperature is 
mapped onto a bosonic system at finite temperature. Equation (|41|) does not determine the 
exact values of the parameters u and T, but knowing the value of the fraction is sufficient 
to fully determine the trace. This can be physically explained by the fact that the original 
spin system is fully randomized at infinite temperature; the number of degrees of freedom 
necessary to fully characterize it is equal to the number of the generators <S&, namely, 3. 
But for the harmonic oscillator, we have only two generators; the z component is found to 
be equivalent to the position variable of another harmonic oscillator that is at zero tem- 
perature. This discrepancy is removed by the value In 3 of the fraction ftw/ksT which is 
clearly the entropy of a fully mixed three-level system. We have applied the formalism to 
the Heisenberg XY model at finite temperature, and found that it is equivalent to a single 
mode harmonic oscillator in thermal equilibrium which can be assigned an effective temper- 
ature that depends logarithmically on the physical one. Moreover, we have determined the 
exact rule for mapping any function of the spin operators to the bosonic counter part. The 
presented results in this paper are a first attempt and further investigation needs to be done 
to gain more insight into other possible applications. 
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